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A Fast Algorithm for Three-Dimensional Layers 
of Maxima Problem 



Yakov Nekrich* 



Abstract. We show that the three-dimensional layers-of-maxima prob- 
lem can be solved in o(n log n) time in the word RAM model. Our al- 



gorithm runs in 0(n(loglogn) ) deterministic time or 0(n(loglogn) ) 
>^r *i expected time and uses 0{n) space. We also describe a deterministic al- 

gorithm that uses optimal 0{n) space and solves the three-dimensional 
layers-of-maxima problem in 0(n log n) time in the pointer machine 
model. 



m 

C/^ ■ 1 Introduction 

ly^ ' A point p dominates a point q if each coordinate of p is larger than or equals to 

O . the corresponding coordinate of q. A point p is a maximum point in a set S if no 

point of S dominates p. The maxima set of S is the set of all maximum points 

psj . in S. In the layers-of-maxima problem we assign points of a set S to layers Si, 

^ ' « > 1, according to the dominance relation: The first layer of S is defined as 

CO , the maxima set of S, the layer 2 of 5* is the maxima set of S' \ 5*1, and the i-th 

layer of S is the maxima set oi S \ (U'"^^^). In this paper we show that the 
three-dimensional layers-of-maxima problem can be solved in o(n log n) time. 

Previous and Related Work. The algorithm of Kung, Luccio, and 
Preparata [53] finds the maxima set of a set S in O(rilogn) time for d = 2 or 

(^ ' d = 3 dimensions and 0(n log ~ n) time for d > A dimensions. The algorithm of 

Gabow, Bentley, and Tarjan [16] finds the maxima set in 0(n log ~ n log log n) 
time for d > 4 dimensions. Very recently, Chan, Larsen, and Patra§cu [TTj de- 
scribed a randomized algorithm that solves the d-dimensional maxima problem 
(i.e., finds the maxima set) for d > 4 in 0{n\og ~^ n) time. Numerous works are 
devoted to variants of the maxima problem in different computational models 
and settings: In [8], the authors describe a solution for the three-dimensional 
maxima problem in the cache-oblivious model. Output-sensitive algorithms and 
algorithms that find the maxima for a random set of points are described 
in [7I13I18I22] . The two-dimensional problem of maintaining the maxima set 
under insertions and deletions is considered in |21j ; the problem of maintaining 
the maxima set for moving points is considered in [151 . 

The general layers-of-maxima problem appears to be more difficult than 
the problem of finding the maxima set. The three-dimensional layers-of-maxima 
problem can be solved in O(nlognloglogn) time [T] using dynamic fractional 
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cascading [511 . The algorithm of Buchsbaum and Goodrich [9^ runs in 0{n log n) 
time and uses 0{n log n log log n) space. Giyora and Kaplan [T7] described a data 
structure for point location in a dynamic set of horizontal segments and showed 
how it can be combined with the approach of 9 to solve the three-dimensional 
layers-of-maxima problem in 0{nlogn) time and 0{n) space. 

The 0(n log n) time is optimal even if we want to find the maxima set in 
two dimensions |23l provided that we work in the infinite-precision computation 
model in which input values, i.e. point coordinates, can be manipulated with 
algebraic operations and compared. On the other hand, it is well known that 
it is possible to achieve o(n log n) time (resp. o(logn) time for searching in a 
data structure) for many one-dimensional as well as for some multi-dimensional 
problems and data structures in other computational models. For instance, the 
grid model, that assumes all coordinates to be integers in the range [1, U] for a 
parameter [/, was extensively studied in computational geometry. Examples of 
problems that can be solved efficiently in the grid model are orthogonal range 
reporting queries |25| and point location queries in a two- and three-dimensional 
rectangular subdivisons [5]. In fact, we can use standard techniques to show 
that these queries can be answered in o(log n) time when all coordinates are ar- 
bitrary integers. Recently, a number of other important geometric problems was 
shown to be solvable in o(nlogn) time (resp. in o(logn) time) in the word RAM 
model. An incomplete lis1o includes Voronoi diagrams and three-dimensional 
convex hulls in 0{n ■ 2'^(^'°siog")) time [T^, two-dimensional point location in 
0(logn/loglogn) time [26110) . and dynamic convex hull in 0(logn/loglogn) 
time [H]. Results for the word RAM model are important because they help us 
better understand the structure and relative complexity of different problems 
and demonstrate how geometric information can be analyzed in algorithmically 
useful ways. 

Our Results. In this paper we show that the three-dimensional layers-of- 
maxima problem can be solved in 0(n(loglogn)'^) deterministic time and 0{n) 
space in the word RAM model. If randomization is allowed, our algorithm runs in 
0(n(loglog?T.)^) expected time. For comparison, the fastest known deterministic 
linear space sorting algorithm runs in O(nloglogn) time [19[ . Our result is 
valid in the word RAM computation model, but the time-consuming operations, 
such as multiplications, are only used during the pre-processing step when we 
sort points by coordinates (see section [5]). For instance, if all points are on the 
n X n X n grid, then our algorithm uses exactly the same model as [25[ or [5]. 

We also describe an algorithm that uses 0{n) space and solves the three- 
dimensional layers-of-maxima problem in optimal O(nlogTi) time in the pointer 
machine model [27]. The result of Giyora and Kaplan [17] that achieved the 
same space and time bounds is valid only in the RAM model. Thus we present 
the first algorithm that solves the three-dimensional layers-of-maxima problem 
in optimal time and space in the pointer machine model. 



^ We note that problems in this list are more difficult than the layers-of-meixima 
problem because in our case we process a set of ajds-parallel segments. 



Overview. Our solution, as well as the previous results, is based on the 
sweep plane algorithm of 9^ described in section [2] The sweep plane algorithm 
assigns points to layers by answering for each p € S a point location query in 
a dynamically maintained staircase subdivision. We observe that general data 
structures for point location in a set of horizontal segments cannot be used to 
obtain an o{n log n) time solution. Even in the word RAM model, no dynamic 
data structure that supports both queries and updates in o(log n) time is known. 
Moreover, by the lower bound of [5] any data structure for a dynamic set of 
horizontal segments needs l2(logn/loglogn) time to answer a point location 
query. We achieve a significantly better result using the methods described below. 

In section [3] we describe the data structure for point location in a staircase 
subdivision that supports queries in 0((loglogn)'^) time and updatef|j in poly- 
logarithmic time per segment. This result may be of interest on its own. 

The data structure of section [3] is not sufficient to obtain the desired runtime 
and space usage mainly due to high costs of update operations. To reduce the 
update time and space usage, we construct auxiliary staircases Bi, such that: 1. 
the total number of segments in Bi and the total number of updates is 0{n/d) 
for a parameter d — log '^^■' n; 2. locating a point p among staircases Bi gives 
us an approximate location of p among the original staircases A4i (up to 0{d) 
staircases). An efficient method for maintaining staircases Bi, described in sec- 
tion m is the most technically challenging part of our construction. In section [5] 
we show how the data structure of section [3] can be combined with the auxiliary 
staircases approach to obtain an 0(T7,(loglogn)^) time algorithm. We also sketch 
how the same approach enables us to obtain an 0{n\ogn) time and 0{n) space 
algorithm in the pointer machine model. 

2 Sweep Plane Algorithm 

Our algorithm is based on the three-dimensional sweep method that is also used 
in [S] . We move the plane parallel to the xy plane |f| from z = +oo to z = and 
maintain the following invariant: when the z-coordinate of the plane equals v all 
points p with p.z > v are assigned to their layers of maxima. Here and further 
p.x,p.y, and p.z denote the x,- y-, and z-coordinates of a point p. Let Si{v) be the 
set of points q that belong to the i-th layer of maxima such that q.z > w; let Pi{v) 
denote the projection of Si{v) on the sweep plane, Pi{v) = {tt{p) \p S Si{v)} 
where 7r(p) denotes the projection of a point p on the xy-plane. For each value of 
V maximal points oi Pi{v) form a staircase A4i; see Fig.[TJ When the z-coordinate 
of the sweep plane is changed from w -|- 1 to u, we assign all points with p.z = v 
to their layers of maxima. If 7r(p), such that p.z = v, is dominated by a point 
from Pi{v + 1), then p belongs to the j'-th layer of maxima and j > i. If 7r(p), 
such that p.z — v, dominates a point on Pk{v -\- 1), then p belongs to the j-th 
layer of maxima and j < k. We observe that 7r(p) dominates Pi{v + 1) if and only 



^ We will describe update operations supported by our data structure in sections [5] 

andO 
^ We assume that all points have positive coordinates. 



if the staircase A4i is dominated by p, i.e., the vertical ray shot from p in —y 
direction passes through Aii. Hence, the point p belongs to the layer i, such that 
Tr{p) is between the staircase M-i-i and the staircase Aii. This means that we 
can assign a point to its layer by answering a point location query in a staircase 
subdivision. When all p with p.z — v are assigned to their layers, staircases are 
updated. 

Thus to solve the layers of maxima problem, we examine points in the de- 
scending order of their z-coordinates. For each v, such that there is at least 
one p with p.z = v, we proceed as follows: for every p with p.z = v operation 
locate(p) identifies the staircase A4i immediately below n{p). If the first stair- 
case below n{p) has index i (7r(p) may also lie on Aii), then p is assigned to the 
i-th layer of maxima; if 7r(p) is below the lowest staircase Aij, then p is assigned 
to the new layer j -\- 1. When all points with p.z — v are assigned to their layers, 
the staircases are updated. All points p such that p.z = v are examined in the 
ascending order of their ^-coordinates. If a point p with p.z = w is assigned 
to layer i, we perform operation replace(p, i) that removes all points of Aii 
dominated by p and inserts p into A4i. If the staircase i does not exist, then 
instead of replace(p, i) we perform the operation new(p, i); new(p, i) creates a 
new staircase A4i that consists of one horizontal segment h with left endpoint 
(0,p.y) and right endpoint 7r(p) and one vertical segment t with upper endpoint 
7r(p) and lower endpoint {0,p.x). See Fig. [T] for an example. 
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Fig. 1. Points a, b, c, d, and e have the same z-coordinate. (a) Points a, b and c 
are assigned to layer 2, d is assigned to layer 3, and e is assigned to a new layer 
6. (b) Staircases after operations replace(a, 2), replace(6, 2), replace(c, 2), 
replace((i, 3), and new(e). Observe that b is not the endpoint of a segment in 
the staircase M2 after updates. 



We can reduce the general layers of maxima problem to the problem in the 
universe of size 0{n) using the reduction to rank space technique [25116) . The 
rank of an element e € S is defined as the number of elements in S that are 
smaller than e: rank(e, 5) — \{a € S\a < e}\; clearly, rank(e, 5) < \S\. For a 



point p = {p.x,p.y,p.z), p € S, let t{p) = {raiiklp.x, Sx) + 1, rank(p.y, 5j^) + 
l,rank(p.z, 5*2) + 1). Let 5' = {t{p) \p G S}. Coordinates of all points in S' 
belong to range [1, n]. A point p dominates a point q if and only if rajik(p.a;, Sx) > 
ra.Tik{q.x, Sx), rajik(p.y, 5j^) > Tank{q.y, Sy), and rank(p.z,5z) > rank(g.z, S'z) 
where Sx, Sy, Sz are sets of x-, y-, and z-coordinates of points in S. Hence if a 
point p' G S" is assigned to the i-th layer of maxima of S' , then T~'^{p') belongs 
to the i-th layer of maxima of S. We can find ranks of x— , y— , and z— coordinates 
of every point by sorting Sx, Sy, and S^- Using the sorting algorithm of flO , S^, 
Sy, and Sz can be sorted in O(nloglogn) time and 0{n) space. Thus the layers 
of maxima problem can be reduced to the special case when all point coordinates 
are bounded by 0{n) in 0(n log log n) time. 



3 Fast Queries, Slow Updates 

In this section we describe a data structure that supports locate(q) in 
0((loglogn)^) time and update operations replace((7, i) and new(g,i) in 
0(logn(loglogn)^) time per segment. We will store horizontal segments of all 
staircases in a data structure that supports ray shooting queries: given a query 
point q identify the first segment s crossed by a vertical ray that is shot from q 
in —y direction; in this case we will say that the segment s precedes q (or s is the 
predecessor segment of q). In the rest of this paper, segments will denote hori- 
zontal segments. Identifying the segment that precedes q is (almost) equivalent 
to answering a query locate(g). Operation replace(g, i) corresponds to a dele- 
tion of all horizontal segments dominated by q and an insertion of at most two 
horizontal segments, see Fig [1] Operation new(g, i) corresponds to an insertion 
of a new segment. 

Our data structure is a binary tree on x-coordinates and segments are stored 
in one-dimensional secondary structures in tree nodes. The main idea of our 
approach is to achieve fast query time by binary search of the root-to-leaf path: 
using properties of staircases, we can determine in 0((loglogn)^) time whether 
the predecessor segment of a point q is stored in the ancestor of a node v or in 
the descendant of a node v for any node v on the path from the root to q.x. Our 
approach is similar to the data structure of [5] , but we need additional techniques 
to support updates. 

For a horizontal segment s, we denote by start(s) and end(s) the x- 
coordinates of its left and right endpoints respectively; we denote by y(s) the 
^/-coordinate of all points of s. An integer e G S* precedes (follows) an integer x 
in S if e is the largest (smallest) element in S, such that e < x {e > x). Let H be 
a set of segments and let Hy be the set of y-coordinates of segments in H. We 
say that s £ H precedes (follows) an integer e if the j/-coordinate of s precedes 
(follows) e in Hy. Thus a segment that precedes a point g is a segment that 
precedes q.y in the set of all segments that intersect the vertical line x — q.x. 

We construct a balanced binary tree T of height log n on the set of all possible 
^-coordinates, i.e., n leaves of T correspond to integers in [1, n]. The range of a 



node V is the interval rng(w) ~ [left(t;),right(u)] where left(w) and right(w) 
are leftmost and rightmost leaf descendants of v. 

We say that a segment s spans a node v if start(s) < lef t(w) < right(w) < 
end(s); a segment r belongs to a node v if left(ti) < start(s) < end(s) < 
right(w). A segment s ^-cuts a node w if s intersects the vertical line x — lef t(i;), 
but s does not span v, i.e., start(s) < left(u) and end(s) < right(t;); a 
segment s r-cuts a node w if s intersects the vertical line x = right(t;) but s does 
not span v, i.e., start(s) > left(i;) and end(s) > right(u). A segment s such 
that [start(s), end(s)] nrng(i;) ^ either cuts v, or spans v, or belongs to v. We 
store y-coordinates of all segments that /-cut (r-cut) a node d in a data structure 
^v (Ti-v)- Using exponential trees 0], we can implement £„ and TZy in linear 
space, so that one-dimensional searching (i.e. predecessor and successor queries) 
is supported in 0((log log n)^) time. Since a segment cuts 0(log n) nodes (at most 
two nodes on each tree level), all £.„ and TZy use 0(?ilogn) space. We denote by 
index(s) the index of the staircase Mi that contains s, i.e., s € A^index(s)- The 
following simple properties are important for the search procedure: 

Fact 1 Suppose that an arbitrary vertical line cuts staircases Aii and A4j, i < j , 
in points p and q respectively. Then p.y > q.y because staircases do not cross. 

Fact 2 For any two points p and q on a staircase Mi, if p.x < q.x, then p.y > 

q-y 

Fact 3 Given a staircase Mi and a point p, we can determine whether Mi is 
below or above p and find the segment s € Mi such that p.x e [start(s), end(s)] 
in 0((loglogn)^) time. The data structure Di that supports such queries uses 
linear space and supports finger updates in 0{1) time. 

Proof: The data structure Di contains x-coordinates of all segment endpoints of 
Mi. Di is implemented as an exponential tree so that it uses 0{n) space. Using 
Di we can identify s G Mi such that p.x G [start(s), end(s)] in 0((loglog7T,)^) 
time; Mi is below p if and only if s is below p. D 

Using Fact [3] we can determine whether a segment s precedes a point q in 
0((loglogn)^) time: Suppose that s belongs to a staircase Mi. Then s is the 
predecessor segment of q iff q.x e [start(s), end(s)], q.y > y(s) and the staircase 
Mi-i is above q. 

We can use these properties and data structures Cy and TZy to determine whether 
a segment b that precedes a point q spans a node v, belongs to a node v, or cuts 
a node v. If the segment b we are looking for spans v, then it cuts an ancestor 
oi v; if that segment belongs to v, then it cuts a descendant of v. Hence, we can 
apply binary search and find in O(loglogn) iterations the node / such that the 
predecessor segment of q cuts /. Observe that in some situations there may be 
no staircase Mi below q, see Fig[3]for an example. To deal with such situations, 
we insert a dummy segment Sd with left endpoint (1, 0) and right endpoint (n, 0); 
we set index(sd) = +oo and store Sd in the data structure Cy^ where vq is the 
root of T. 



Let Ix be the leaf in which the predecessor of q.x is stored. We wih use 
variables I, u and v to guide the search for the node /. Initially we set / = l^ 
and u is the root of T- We set v to be the middle node between u and /: if the 
path between u and I consists of h edges, then the path from u to v consists of 
\h/2\ edges and v is an ancestor of I. 

Let r and s denote the segments in £„ that precede and follow q.y. If there 
is no segment s in C^, with y{s) > q.y, then we set s = NULL. If there is no 
segment r in C^ with y{r) < q.y, then we set r — NULL. We can find both r 
and s in 0((loglogn)^) time. If the segment r ^ NULL, we check whether the 
staircase A^index(r) contains the predecessor segment of q; by Fact [31 this can 
be done in 0((loglogn)^) time. If 7Windex(r) contains the predecessor segment 
of q, the search is completed. Otherwise, the staircase -Mindex(r)-i is below q or 
r = NULL. In this case we find the segment r' that precedes q.y in 7?,^. If r' is 
not the predecessor segment of q or r' = NULL, then the predecessor segment 
of q either spans v or belongs to v. We distinguish between the following two 
cases: 

1. The segment s ^ NULL and the staircase that contains s is below q. By 
Fact [H a vertical line x = q.x will cross the staircase of s before it will cross a 
staircase Mi, i > index(s). Hence, a segment that spans v and belongs to the 
staircase Mi, i > index(s), cannot be the predecessor segment of q. If a segment 
t spans V and index(i) < index(s), then the y-coordinate of t is larger than the 
2/-coordinate of s by Fact[TJ Since y{t) > y{s) and y{s) > q.y, the segment t is 
above q. Thus no segment that spans v can be the predecessor of q. 

2. The staircase that contains s is above q ot s = NULL. If r exists, the 
staircase A^imiox(r)-i is below q. Hence, the predecessor segment of q belongs 
to a staircascQ Mi, index(s) < i < index(r) — 1. Since each staircase Mi, 
index(s) < i < index(r) — 1, contains a segment that spans v, the predecessor 
segment of s is a segment that spans v. If r does not exist, then every segment 
below the point q spans the node v. Hence, the predecessor segment of s spans 
V. See Fig. [5] for an example. 

If the predecessor segment spans v, we search for / among ancestors of v; if 
the predecessor segment belongs to v, we search for / among descendants of v. 
Hence, we set I = v in case 2, and we set u = u in case 1. Then, we set v to be 
the middle node between u and / and examine the new node v. Since we examine 
O(loglogn) nodes and spend 0((loglogri)^) time in each node, the total query 
time is 0((loglogn)'^). 

If the predecessor segment is the dummy segment Sd, then there is no hori- 
zontal segment of any Mi below q. In this case we must identify the staircase to 
the left of q.x. Let rrii denote the rightmost point on the staircase Mi, i.e., rrii 
is a point on Mi such that rrii.y — 0. Then q is between staircases Mi-i and 
Mi, such that rrii.x < q.x < irii^i.x. We can find rrii in 0((loglogn)^) time. 

When a segment s is deleted, we delete it from the corresponding data struc- 
ture Di. We also delete s from all data structures C^ and 7?,^ for all nodes v and 
w, such that s /-cuts v (respectively r-cuts w). Since a segment cuts O(logn) 



To simplify the description, we assume that index(s) = if s = NULL. 
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Fig. 2. Search procedure in a node v. Staircases are denoted by their indexes. 
Figures (a) and (b) correspond to cases 1 and 2 respectively. The case when 
the predecessor segment belongs to A^index(r) is shown on Fig. (c). 



nodes and exponential trees support updates in O((loglog7i)^) time, a deletion 
takes 0(log n(log log n)^) time. Insertions are supported in the same wajO. Oper- 
ation new(g, I) is implemented by inserting a segment with endpoints (0, q.y) and 
{q.x,q.y) into T, incrementing by one the number of staircases I, and creating 
the new data structure Di. To implement replace(g,i) we delete the segments 
"covered" by q from T and Di and insert the new segment (or two new segments) 
into T and Di. 

Lemma 1. We can store n horizontal staircase segments with endpoints on 
n X n grid in a 0{n\ogn) space data structure that answers ray shoot- 
ing queries in 0((loglogn)'^) time and supports operation replace(q, i) in 
0(TOlogn(loglogn)^) time where m is the number of segments inserted into 
and deleted from the staircase A4i, and operation new(g) in O(logn(loglogn)^) 
time. 

The data structure of Lemma [T] is deterministic. We can further improve the 
query time if randomization is allowed. 

Fact 4 Given a staircase M.i and a point p, we can determine whether A4i is 
below or above p and find the segment s £ A4i such that p. x G [start(s), end(s)] 
in 0((loglogn)) time. The data structure Di that supports such queries uses 
linear space and supports finger updates in 0(1) expected time. 

Proof : The data structure is the same as in the proof of Fact [H but we use the 
y-fast tree data structure [28] instead of the exponential tree. D 

Lemma 2. We can store n horizontal staircase segments with endpoints onnxn 
grid in a 0{nlogn) space data structure that answers ray shooting queries in 
0((loglogn)^) time and supports operation Tepl&ce{q,i) in 0{m log nloglogn) 
expected time where m is the number of segments inserted into and deleted from 
the staircase A4i, and operation ne'w{q) in 0(log nloglogn) expected time. 



^ The update time can be slightly improved using fractional cascading and similar 
techniques, but this is not necessary for our presentation. 



Proof: Our data structure is the same as in the proof of Lemma [TJ But we 
implement Di using Fact ID Data structures C^ and TZy are implemented using 
the y-fast tree [2B]- Hence, the search procedure spends 0(loglog7i) time in each 
node of T and a query is answered in 0((loglogn)^) time. D 

Although this is not necessary for further presentation, we can prove a similar 
result for the case when all segment endpoints are on a U x U grid; the query 
time is 0{\og\ogU + (log log n)^) and the update time is 0(log n(loglogn)^) 
per segment. See Appendix D for a proof of this result. 



4 Additional Staircases 

The algorithm in the previous section needs 0(71 log n(log log n)^) time to con- 
struct the layers of maxima; n ray shooting queries can be performed in 
0(n(loglogn)^) time, but 0{n) update operations take 0(nlogn(loglogn)^) 
time. To speed-up the algorithm and improve the space usage, we reduce the 
number of updates and the number of segments in the data structure of Lemma[T] 
to 0(n/log^n). 

Let V denote the data structure of Lemma [1] We construct and maintain a 
new sequence of staircases 81,82,- ■ ■ ,8m, where m < n/d and the parameter 
d will be specified later. All horizontal segments of Si, ... , Bm are stored in V. 
The new staircases satisfy the following conditions: 

1. There are 0(§) horizontal segments in all staircases 8i 

2. X* is updated O(^) times during the execution of the sweep plane algorithm. 

3. For any point q and for any i, if q is between Si_i and Bi, then q is situated 
between Mk and Mk+i for {i - 3/2)d < k < {i + l/2)d. 

Conditions 1 and 2 imply that the data structure V uses 0{n) space and all 
updates of V take 0{n) time if d > logn(loglogn)^. Condition 3 means that 
we can use staircases 8i to guide the search among A4k'- we first identify the 
index i, such that the query point q is between B^+i and 8i, and then locate 
q in A^(i-3/2)dj • ■ • : -'^(i+i/2)d- It is not difficult to construct Bi that satisfy 
conditions 1 and 3. The challenging part is maintaining the staircases Bi with a 
small number of updates. 

Lemma 3. The total number of inserted and deleted segments in all 8i is 0(4). 
The number of segments stored in Bi is O(^). 

We describe how staircases can be maintained and prove Lemma [3] in Appendix 
B. 



5 Efficient Algorithms for the Layers-of-Maxima Problem 

Word RAM ModeL To conclude the description of our main algorithm, we 
need the following simple 



Lemma 4. Using a 0{m) space data structure, we can locate a point in a group 
of d staircases A4j, A4j+i, . . . ,Aij+d in 0{\ogd ■ (log log m)^) time, where m is 
the number of segments in A4j, A^j+i, . . . , Aij+d- An operation replace(g, i) is 
supported in 0((loglog?Ti)^ +rnq) time, where ruq is the number of inserted and 
deleted segments in the staircase M.i, j < i < i + d. 

Proof: We can use Fact [3] to determine whether a staircase is above or below a 
staircase Mk for any j < k < j + d. Hence, we can locate a point in 0{\ogd ■ 
(log log m)^) time by a binary search among d staircases. D 

We set d = log n. The data structure Fi contains all segments of staircases 
M(^i-i)d+i, M[i-i)d+2, ■ ■ ■ , Mtd for i = l,2,...,j, where j = [l/d\ and I is the 
highest index of a staircase; the data structure -Fj+i contains all segments of 
staircases Mjd+i, ■ ■ ■ ,M.i. We can locate a point q in each Fi in 0((loglogn)^) 
time by Lemma HI Since each staircase belongs to one data structure, all Fi use 
0{n) space. We also maintain additional staircases Bi as described in section 21 
All segments of all staircases Bi are stored in the data structure V of Lemma [TJ 
since T) contains 0{n/d) segments, the space usage of V is 0{n). 

Now we can describe how operations locate, replace, new can be imple- 
mented in 0((loglogn)^) time per segment. 

— locate((7): We find the index fc, such that q is between Bk-i and Bk in 
0((loglogn)^) time. As described in section IH q is between Mkd+g and 
M{k-i)d-g- Hence, we can use data structures Fk+i, Fk, and Fk-i to identify 
j such that q is between A4j and Aij+i. Searching Fk+i, Fk, and Fk-i takes 
0((loglogn)^) time, and the total time for locate((j) is O((loglogn)^). 

— replace(g, i): let niq be the number of inserted and deleted segments. The 
data structure F^i/dj can be updated in 0{mq + (log log n)^) time. We may 
also have to update B^i^d] ■, ^li/d\+i: ^nd the data structure V. 

— new(g, Z): If Z = kd + 1 for some fc, a new data structure Fk^i is created. 
We add the horizontal segment of the new staircase into the data structure 
Ffe_|_i. If / — kd, we create a new staircase Bk and add the segments of Bk 
into the data structure T). 

There are 0{n/d) update operations on the data structure T) that can be per- 
formed in 0((ri/(i) log n(log log n)'^) — 0{n) time. If we ignore the time to 
update T>, then replace(g,i) takes 0(?7i^(loglogn)^) time and new(q, /) takes 
0((loglogn)^) time. Since X^„g5?7ig = 0{n) and new(g, Z) is performed at most 
n times, the algorithm runs in 0(ri(loglogn)^) time. We thus obtain the main 
result of this paper. 

Theorem 1. The three-dimensional layers- of -maxima problem can be solved in 
0(n(loglogn)^) deterministic time in the word RAM model. The space usage of 
the algorithm is 0{n). 

If we use Fact U instead of Fact |3] in the proof of Lemma 2] and Lemma [3] instead 
of Lemma [1] in the proof of Theorem [ll we obtain a slightly better randomized 
algorithm. 



Theorem 2. The three-dimensional layers- of -maxima problem can be solved in 
0(n(loglogn)^) expected time. The space usage of the algorithm is Oin). 

Pointer Machine Model. We can apply the idea of additional staircases 
to obtain an 0(n log n) algorithm in the pointer machine model. This time, we 
set d = log n and maintain additional staircases Bi as described in section SI 
Horizontal segments of all Bi are stored in the data structure T> of Giyora and 
Kaplan [17^ that uses 0(?Tilog'^ ?n) space and supports queries and updates in 
0(logr7T.) and 0(log ^'^ to) time respectively, where to is the number of segments 
in all Bi and e is an arbitrarily small positive constant. Using dynamic fractional 
cascading [M], we can implement Fi so that Fi uses linear space and answers 
queries in 0(logn + log log n log d) = O(logn) time. Updates are supported in 
0(log n) time; details will be given in the full version of this paper. Using 2? and 
Fi, we can implement the sweep plane algorithm in the same way as described 
in the first part of this section. The space usage of all data structures Fi is 0(n), 
and all updates of Fi take O(nlogn) time. By Lemma [31 the data structure V 
is updated 0(n/ log n) times; hence all updates of V take O(nlogn) time. The 
space usage of 2? is 0{m\og^ m) = 0(n). Each new point is located by answering 
one query to V and at most three queries to Fi; hence, a new point is assigned 
to its layer of maxima in O(logn) time. 

Theorem 3. A three-dimensional layers- of-maxima problem can be solved in 
0{n\ogn) time in the pointer machine model. The space usage of the algorithm 
is 0{n). 
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Appendix A. Figures 
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Fig. 3. Example of a just constructed additional staircase Bi for d — 6. The 
staircase Bi is shown with dashed red lines. Staircases are denoted by their 
indexes. 



^ 



Fig. 4. There are no staircases below p and q. 
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Fig. 5. Maintaining Invariant |4l For simplicity only several segments of Aiid, 
■Mid-g, and Bi are shown; Bi is drawn in red color, (a) & (b): Two empty 
segments are replaced with one (possibly empty) segment t' . (c) & (d): Two 
empty segments are replaced with one non-empty segment i' and one empty 
segment t" . 



Appendix B. Proof of Lemma [3] 



In the first part of this section we describe the construction procedure of a 
boundary Bi. Then, we will prove some facts about Bi and describe the update 
procedure. In the last part of this section we will prove that all Bi are updated 
0(1) times for d updates of A^;. 

Construction of Additional Staircases. We construct one staircase Bi 
for d staircases A^(i_i)d+i, • ■ • , Mid- Let p be the starting point of the staircase 
Mid, i-e., p & Mi and p.x = 1. The staircase Bi is the path traced by p as we 
alternatively move p in the -l-a; and —y direction until it hits the cc-axis. 

A segment s covers a point p if start(2;) < p.x < end(s). A segment r is 
related to a segment s if s covers the left endpoint of r; a segment s covers a 
segment r if start(s) < start(r) and end(r) < end(s). A point p dominates a 
segment s Up dominates the left endpoint of s. A segment s follows the segment 
r in a staircase Bi or Mi (resp. r precedes s) if both r and s belong to the same 
staircase and end(?') = start(s). 



Let g — d/2. For convenience we assume that each point q ^ S has even x- 
coordinate. This is achieved by replacing each point q — {q.x, q.y) with a point 
q' = {2q.x, q.y). Endpoints of ah segments of Bi will have odd x-coordinates. The 
set Gi contains all segments of A4id-g, ■ ■ ■ , Mid- The staircase is constructed by 
repeating the following steps until p hits the a;-axis or the x-coordinate of p is 
maximal possible, i.e. until p.y = or p.x = 2n: 

(1) We movep in the +x direction untilp cuts A4id~-g, i-e until p.x = start(s) + l 
for a segment s G A4id-g such that y(s) < p.y 

(2) If p.x < 2n, we move p in —y direction until it hits a segment of A4id or 
p.y = 0. 

Observe that at the beginning of step (1) the point p always belongs to a horizon- 
tal segment of Aiid- Hence, a point on A4id+i does not dominate a segment of Bi. 
Since each horizontal segment of Bi cuts Mid^g it also cuts Mid-j, < j < g. 
Hence, there are at least g segments of Gi related to each horizontal segment of 
Bi and the total number of segments in all Bi is 0{-). An example of a (just 
constructed) additional staircase is shown on Fig. [H 

Updates. When we update a staircase Aiid+j for g/2 > j > —5/2 by operation 
replace, the staircase is moved in the north-east direction. As a result, a point 
on a staircase A4id+j, j > 0, may dominate a segment of Bi. Therefore we main- 
tain a weaker property: no segment of A4id+g dominates Bi and each point of 
Bi is dominated by a point on A4id~g- Our goal is to update Bi 0(1) times for 
f^{g) updates of Mj (in average). We achieve this by maintaining the following 
invariants 

Invariant 1 Each segment s e Bi is dominated by the right endpoint of a seg- 
ment r G Mid- 

Invariant 2 No point of Bi is dominated by a point of Aiid+g/2- 

Invariant 3 No segment s G Bi cuts A4id-g+i - 

We say that a segment s is empty if it does not cut Mid- g/2- 

Invariant 4 If a segment S2 follows si in Bi, then either S2 or si is not empty. 

If Invariants [I] and [3] are true when Bi is constructed, they will not be violated 
after updates oiMid,---,Mid-g- We update Bi if Invariants [2] or |4] are violated: 
If a segment s G Bi, such that s was not empty when s was inserted into Bi , does 
not cut Mid-g/2 after an operation replace((7, id — g/2), we call the procedure 
Rectify(Si,s) that will be described later in this section. If a segment s of Bi 
is dominated by a point of Mid+g/2 after replace((7, id + g/2), we also call the 
procedure Rectif y(i3i, s). 

Fact b If a point q dominates a segment of Mid- j, then q dominates at least 
one segment of Mid- k for each k < j. 

Fact 6 If a point q dominates more than two segments of Bi, then q dominates 
at least one segment of Mid-j for each j — 1,2, . . . ,g/2. 



Proof: If q dominates three segments of Bi then q dominates the right endpoint 
of at least one non-empty segment s € Bi. Since s cuts Aiid~g/2i the right 
endpoint of s dominates a segment of A^id-g/2- Hence, the right endpoint of s 
dominates at least one segment of Mid-j for j = 1, 2, . . . ,5/2 — 1 by Fact [S) 
Since q dominates the right endpoint of s, g also dominates at least one segment 
of Mid-j for each j = 1, 2, . . . , g/2. D 

Fact 7 Any point q on Mid+j, j > 0, dominates at most two segments of Bi. 

Proof: Suppose that a point q on M-id+j dominates more than two segments 
of Bi. Then, there is a point q' on Mid that also dominates more than two 
segments of Bi. By Fact [SI q' dominates a segment of Mid-j for each j = 
1, 2, . . . , 5/2. Since a point on Mid cannot dominate a point on Mid-g/21 we 
obtain a contradiction. D 

Fact [3 which is a corollary of Invariant |31 guarantees us that each operation 
replace((7, id + j) such that q dominates Bi affects at most two segments of Bi. 
This will be important in our analysis of the number of updates of Bi. Now we 
are ready to describe the update procedure. 

The procedure Rectif y(Si, s) deletes a segment s and a number of preceding 
and following segments and replaces them with new segments. We say that a 
segment s' is the child of s if s was removed by an operation replace((7), such 
that q is the right endpoint of s'; s' is a descendant of s if s' is a child of s 
or a descendant of a child of s. Let sq be the segment that precedes s in Bi. 
Let si, S2, . . . be segments of Bi such that si follows s and Si follows Si-i for 
i > 1. Suppose that s contained the right endpoint of a segment r^ ^ Mid that 
belonged to Mid when s was inserted into Bi, and let r„ be the descendant of 
Tb that belongs to Mid when the procedure Rectify(Bi, s) is performed. By 
Fact [71 r„ may dominate the segment sq that precedes s in 6^ , but r„ does not 
dominate the segment that precedes sq in Bi. Let rmax be a descendant of a 
segment r G Alid related to s with the largest x-coordinate of its right endpoint. 
By Fact [3 rmax may dominate si and S2 but it cannot dominate S3. 

Now we must decide which segments are to be deleted from Bi and how to 
construct new segments. We delete segments that are dominated by r„ or Tmax- 
As shown above, there are at most three such segments (except of s itself). If s/, 
/ < 2, is the last segment dominated by rmax, we may also remove some segments 
that follow Sf. But our guarantee is that all removed segments s/+i, ■ . ■ Sm do 
not cut Mid-^g/A- We insert new segments into Bi by moving a point p in -l-a; 
and —y directions. A more detailed description follows. 

To simplify the description, we will use set Vi that contains some horizontal 
segments that currently belong to Mid and some segments that belonged to 
Mid but are already deleted. When a staircase Bi is constructed, Vi contains 
all horizontal segments of Mid- When the procedure Rectify(;Bi, s) is called, 
we delete all segments of Vi dominated by r„ or Tmax and insert all segments 
r e Mid, such that end(rK) < start(r) < start(r,„ax)- Segments of Vi are used 
to "bound the staircase Bi from below" , i.e. the left endpoint of each horizontal 



segment in Bi belongs to a segment from. Vi- 

(1) Let p be the point on Bi such that p.y = y(Tv)- This is the left endpoint of 
the first inserted segment of Rectif y(Bi, s). 

(2) We move p in the +x direction until p.x ~ end(rmax) or p cuts M.id~g- While 
p.x < end(r„iax), we repeat the following steps: we move p in the —y direction 
until it hits "new" Mid', then, we move p in +x direction until it cuts A4id-g- 
Observe that all horizontal segments inserted in step 2 cut Aiid-g- 

(3) When p.x = end(rniax)5 we move p in +x direction until it cuts A4id-3g/4 
and p.x > end(r,iiax) + 1- Suppose that now start(sm) < p.x < end(sm) for 
some Sm in Bi. We continue moving p in +x direction until p.x — end(sm) or p 
cuts A4id-g- Then, we move p in —y direction until p hits Vi 

(4) We insert a new segment t instead of the segment Sm- It is possible that now 
p.y < y{sm)- We move p in +x direction until p.x = end(s„i) and move p in —y 
direction until p hits Vi. 

(5) Now we must pay attention that Invariant 2] is maintained: all inserted seg- 
ments except of may be the last one are not empty. Let t denote the last in- 
serted segment, and suppose that both t and Sm+i are empty. We can replace 
two empty segments either with one empty segment, or one non-empty and one 
empty segment as follows. We replace t with a new segment t' : the left endpoint 
of t' coincides with the left endpoint of t and either end(i') = end(sm+i) or t' 
cuts A4id-g- If t' cuts Aiid-g, wc rcplacc Sm+i with a new segment t" such that 
start(i") = end(t'), end(i") — end(s,„+i), and y(t") = y{srn+i)- See Fig. [5]for 
an example of step (5). 

Observe that all but one non-empty segments constructed by Rectif y(;Bi, s) cut 
Mid-g- The only exception is the segment constructed in step (3) that cuts 
Mid-3g/4- We prove in Appendix C that the total number of updates of Bi 
during the execution of the sweep plane algorithm is O(^). 
This completes the Proof of Lemma 121 

Appendix C. Analysis of Update Operations for Additional 
Staircases 

We will show below that the data structure V is updated O(-) times during the 
execution of the sweep-plane algorithm. First, we will estimate the number of 
deleted segments. We will estimate the number of insertions in the end of this 
section. We assign c credit points to each segment of Mid-j and 3c credit points 
to every segment of Aiid+j for j = 1,2, ... ,g and c = 24. Insertion of a new 
segment into Bi is free and deletion costs g credit points. 

Every time when we perform operation replace((7,icf -f- j) for j > we 
distribute the credit points of the newly inserted segment r with right endpoint 
q among several segments of Bi . We evenly distribute credits of r among segments 
Sj G B such that either q dominates Sj or q.x > start(sj) and q.y > j/(sj+i) 
where Sj+i denotes the segment that follows Sj in Bi. By Fact [B] there are at 
most three such segments Sj; hence each Sj obtains at least c credits. When a 
segment r of Mid-j, j > 0, is deleted, we assign credits of r to s G Bi, such 



that r is related to s. We say that a segment s is initiahy non-empty if s cuts 
Mid-g/2 when 3 is inserted into Bi. We will show below that we can pay g credit 
points for each deleted segment of Bi and maintain the following property. 

Property 1. Every initially non-empty segment s ^ Bi that does not cut 
M.id-g+k, k > g/A, accumulated at least c ■ k credit points. 
Every segment s € Bi that is dominated by a point on AAid+j accumulated j ■ c 
credit points. 

Proof: Property [T] is obviously true for a just constructed staircase Bi. Suppose 
that Property[T]is true after the procedure Rectif y(Si, s) was called for segments 
of Si / > times. We will show that this property is maintained after the (/-|-1)- 
th call of the procedure Rectify(;Bi, s). If an initially non-empty segment s does 
not cut Mid-g+k' after the /-th call of Rectify is completed, then s accumulated 
k'c credits. If the segment s does not cut Mid-g+k at some point after the /-th 
call of Rectify, then at least k — k' segments of Mid-g+j, k' < j < k, that are 
related to s are already deleted. Hence, s accumulated at least ck' + c{k~k') = ck 
credits. Therefore if an initially non-empty segment s does not cut A4id-g/2i then 
s accumulated cg/2 credits. 

If a point of Aiid+k' dominates a segment s ^ Bi when the /-th call of the 
procedure Rectify is completed, then s has k' ■ c credit points. If a point of 
A4id+k dominates a segment s G Bi, then we performed at least one operation 
replace(q, id + j) such that q dominates s for each k' < j < k. Hence, s 
accumulated fee credits. Therefore if a segment of Bi is dominated by a point on 
■Mid+g/2, then s has cg/2 credits. 

Hence, when we start the procedure Rectif y(Si, s), the segment s has cg/2 
credit points. In addition to s, we may have to remove segments sq,si, S2 because 
a descendant of some segment r £ Mid, such that r was related to s when s 
was constructed, dominates sq, si, or S2- If segments S3, S4, . . . s„i are removed 
by Rectif y(i3i, s), then each non-empty segment among S3, . . . ,Sm does not cut 
Mid-3g/4:- By Property [U every such segment has eg /A credits. Since there are 
at least (m — 4)/2 non-empty segments among S3, ... , Sm, we can use cg{m — 4)/4 
credits accumulated by non-empty segments to remove S4, . . . , Sm- We use eg /A = 
6g credits accumulated by s to remove s and to remove so,si,S2,S3, and Sm-i-i if 
necessary. If a segment s' inserted after the procedure Rectify cuts only Sg/A 
staircases Mid+j, then we transfer to s' the remaining credit points accumulated 
by s. Recall that there is at most one such segment s' that may be inserted during 
step (3) of the update procedure. Since s accumulated at least eg/2 credit points 
and at most cg/i are spent for removing the segments s, sq, si, and S2, the 
segment s' obtains at least cg/4 credit points after the update procedure. 

We must also take care of the segment t resp. segments t' and t" . Since it is 
possible that y{t) < y(sm), t can be dominated by a point of Mid+j for some j > 
when it is constructed. Remaining credit points of segment Sm are transferred 
to t (resp. to t'); if t" is constructed, then credit points of Sm+i are transferred 
to t" . If t' is constructed but t" is not constructed (i.e. if t' replaces both Sm 
and Sm+i), then credits of Sm+i are also transferred to t' . li t is dominated 



by a point of Mid+k for some < k < g/2, then we performed an operation 
replace((7, zrf + j) for each j < k, such that q.y > y{t) and q.x > start(i). 
Since start(i) > start(s,„) and y{t) > y{sm+i), s,„ was assigned c credits for 
each 1 < j < fc. Hence, if t is dominated by A4id+j, then t has at least cj credit 
points. The same is also true for t' and t" . D 

We can conclude from Property [T] that we can always pay g credit points for 
a deleted segment of Bf, hence, the total number of deleted segments is 0{-). 

Let Nf be the number of segments in all Bi when the algorithm is fin- 
ished. Since for every second segment s in Bi there are at least g/2 segments of 
Mid, ■ ■ ■ ,Mid-g related to s, the total number of segments in Bi is 0{ni/g) 
where n^ denotes the total number of segments in Mid,- ■ ■ T-M-id-g- Hence, 
Nf = 0{Y,nJg) = 0{n/g). Clearly Ni = Nf + Nd where Ni is the num- 
ber of inserted segments and Nd is the number of deleted segments. Hence, 
Ni = 0{n/g) and the total number of inserted and deleted segments in all Bi is 
N, + Nd = Oin/g) = Oin/d). 



Appendix D. Staircases on U X U Grid 

Lemma 5. We can store n horizontal staircase segments with endpoints on 
U X U grid in a 0{n\ogn) space data structure that answers ray shooting 
queries in 0{log\ogU + (log log n)^) time and supports operation replace(g,i) 
in O (to log n(loglogn)^) time where m is the number of segments inserted into 
and deleted from the staircase Aii, and operation new(g) in 0(log n(loglogn)^) 
time. 

Proof: Instead of storing point coordinates of segment endpoints in the data 
structure, we store labels of point coordinates: each x— and y-coordinate is 
assigned an x-label (y-label), so that the x-label (y-label) of q is smaller than 
the x- label (y-label) of p if and only if q.x < p.x {q.y < p.y). All labels belong 
to range [l,0(n)] and are maintained using the technique of |20|29) . When a 
new segment is inserted or deleted, 0(log^ n) labels may change, and we have 
to delete and re-insert into data structures those segments whose labels are 
changed. Since each segment is stored in 0(log n) secondary data structures, a 
deleted/inserted segment leads to O(log^n) updates in C^ and TZ^. Hence the 
update time is 0(log'^ ri(loglogri)^). The query procedure is exactly the same as 
in the proof of Lemma [TJ D 



